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Subject:  Time-Frequency  Hop  Coding  Based  Upon 

the  Theory  of  Linear  Congruences 


ABSTRACT 

Time- frequency  hop  codes  are  developed  based  upon  the  theory  of 
linear  congruences.  These  codes  can  be  used  for  multiuser  radar  and 
asynchronous  spread-spectrum  communications  systems.  A  uniform  upper 
bound  is  placed  on  the  crosscorrelation  function  between  any  two  elements 
of  the  code  set.  The  upper  bound  is  minimized  by  choice  of  BT-product 
and  is  shown  to  diminish  as  2/N,  where  N  is  the  number  of  elements  in 
the  code  set.  The  size  and  position  of  spurious  peaks  in  the  auto¬ 
correlation  functions  are  discussed.  The  results  are  extended  to 
narrowband  ambiguity  functions. 
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INTRODUCTION 

In  chis  paper,  we  investigate  a  sequence  of  time-frequency  hop  codes 
which  can  be  used  in  the  areas  of  coherent  multiuser  radar  and  asynchronous 
spread  spectrum  communications.  Unlike  synchronous  communication  systems 
in  which  one  usually  searches  for  code  sets  which  have  mutually  ortho¬ 
gonal  elements, ^  these  applications  require  that  the  entire  crosscorrelation 
function  between  any  two  elements  of  the  set  be  small.  This  is  true  since 
the  exact  time  of  arrival  of  the  received  signal  is  unknown.  Thus,  the 
output  of  the  matched  filter  for  the  correct  signal  must  be  small  when 
any,  or  perhaps  all,  of  the  other  codes  are  present  within  the  received 
signal. 

In  addition  to  the  bound  that  must  be  placed  on  the  crosscorrelation 
functions  in  these  applications,  several  other  constraints  may  also  be 
placed  upon  the  code  set.  If  time-of -arrival  (i.e.,  range  measurement  or 
resolution  of  closely  spaced  targets)  is  of  importance,  then  the  auto¬ 
correlation  functions  for  the  code  words  must  have  pulse  compression 
(2) 

properties.  Thus,  they  must  possess  a  narrow  main  lobe  and  small 
spurious  peaks  away  from  the  main  lobe.  If  the  transmitter,  target  or 
receiver  is  in  motion,  one  must  then  deal  with  Doppler  effects  and,  hence, 
the  associated  ambiguity  functions  for  the  signal  set.  Thus,  one  would 
like  the  uniform  bound  placed  on  the  crosscorrelation  functions  for  the 
code  set  to  extend  to  the  entire  crossambiguity  function  surface.  Finally, 
depending  upon  whether  or  not  these  relative  motions  are  to  be  measured  or 
Ignored  determines  additional  constraints  on  the  autoambiguity  function  for 
the  code  set.  Thus,  if  velocity  measurement  is  to  be  performed,  then  the 
ideal  "thumb  tack"  ambiguity  function  would  be  desired.  However,  should 
relative  velocity  be  unimportant,  then  the  codes  should  possess  Doppler 

(4) 

tolerant  properties  (i.e.,  a  "razor  blade"  autoambiguity  function  with 
minimum  range-Doppler  coupling). 

The  code  set  presented  here,  whose  structure  is  based  upon  the  theory 
of  linear  congruences,  possesses  most  of  the  desirable  properties  described 
above.  Using  the  number-theoretic  properties  of  linear  congruences,  we 
provide  an  upper  bound  on  the  crosscorrelation  function  for  any  two  elements 
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of  che  code  set.  Further,  we  find  an  expression  relating  the  time-bandwidth 
product  occupied  by  the  code  set  to  the  nuaiber  of  code  words  within  the  set 
which  minimizes  the  crosscorrelation  upper  bound.  Here  bandwidth  is  defined 
as  the  range  of  sinusoidal  frequencies  spanned  by  the  set,  ignoring  the 
spectral  skirts.  Since  all  elements  of  the  code  set  span  all  available 
frequencies,  each  code  word  possesses  maximum  pulse  compression  capabilities. 

A  bound  is  also  placed  on  the  amplitude  of  the  spurious  peaks  in  the  auto¬ 
correlation  function  showing  the  trade-offs  involved  in  the  choice  of  the 
time-bandwidth  product,  the  code  set  size,  and  the  amplitude  of  the 
spurious  peaks.  In  addition,  the  positions  of  the  major  spurious  peaks 
are  also  predictable. 

Finally,  since  che  code  words  have  virtually  identical  time  delay  and 
frequency  shift  properties,  these  results  extend  directly  to  the  narrow- 
band  crossambiguity  functions  for  the  code  set. 

II.  THE  CODE  WORD  SETS 

We  consider  a  rectangular  pulse  of  length  T  seconds,  divided  Into  N 
equal  segments  of  length  T/N  seconds.  For  the  sequel,  N  will  be  restricted 
to  be  a  prime  number.  We  place  in  each  segment  of  the  pulse  a  cosine  wave 
whose  frequency  is  one  of  the  N  frequencies  defined  by 

wk  -  “,0  +  ^  ’  k  -  0,  1 . N  -  1  (1) 

where  B  is  defined  as  the  approximate  bandwidth  of  the  signal.  Thus,  each 
code  word  occupies  a  time-bandwidth  product  of  approximately  2  BT.  We 
shall  assume  w »  to  be  sufficiently  high  to  assure  we  are  dealing  with 
analytic  signals.  We  place  one,  and  only  one  frequency,  in  each  time  slot 
by  using  a  linear  congruent ial  relationship  as  follows:  consider  the 
linear  congruence 


k(m,£)  -  ml. (Mod  N) ,  0  <  m,  l  <  N 


(2) 
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III. 


where  k  ia  the  frequency,  ■  ie  the  code  word  Index,  t  the  time  slot,  and 
N  la  a  prime  nuaber. 

As  an  example,  conaider  the  coda  array  generated  for  N  -  5.  This  is 

ahown  In  Figure  1.  Note  that  except  for  the  left  aoat  column,  £  ■  0,  the 

array  la  coapletely  occupied.  The  nuaber  inalde  each  block  la  the  code 

word  Index.  The  code  word  Indexed  by  a  ■  0  la  a  CW  pulae  and  will  be 

dropped  froa  the  aet.  We  obaerve  that  the  code  word  a  ■  1  la  the  same 

(2) 

pulse  Klauder,  et  al.  uaed  to  heur iatlcally  derive  the  chirp  radar  pulse. 
Thus,  this  procedure  ia  a  generalization  of  linear  frequency  modulation. 


k,  frequency  slot 
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Figure  1.  Time-Frequency  Hop  Codes  for  N  -  5 


In  the  sequel,  we  shall  first  review  certain  properties  of  linear 
congruences  in  the  form  of  Equation  (2)  and  then,  using  these  properties, 
establish  the  crosscorrelation  upper  bounds. 

PROPERTIES  OF  LINEAR  CONGRUENCES 

In  this  section,  we  review  the  properties  of  linear  congruences. 

(5) 

For  a  complete  treatment  of  the  theory  of  linear  congruences,  see  Griffen 
from  which  most  of  the  results  presented  here  are  taken. 

Two  Integers  are  congruent  (Mod  N)  if,  and  only  if,  their  difference 
la  dlvisable  by  N  i  0.  The  totality  of  integers  congruent  to  a  given 
integer  for  modulus  N  constitutes  a  residue  class  (Mod  N).  Any  set  of 
) N |  integers  selected  so  that  no  two  of  them  belong  to  the  same  residue 
class  (Mod  N)  foras  a  complete  residue  aystea  (Mod  N).  Finally,  we  define 
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Che  complete  residue  system  (Mod  N)  which  contains  the  elements  of  each 
residue  class  with  smallest  absolute  magnitude  as  the  minimally  complete 
residue  system  (Mod  N).  As  shorthand  notation,  we  use  RC,  CRS,  and  MCRS, 
respectively. 

As  an  example,  for  N  *  5,  we  have 

a)  1,  6,  and  -4  are  elements  of  the  same  RC 

b)  {5,1, 7, -2, 4}  is  a  CRS 

and 

c)  {-2, -1,0, 1,2}  is  the  MCRS. 

In  general,  for  N  prime  and  greater  than  2,  the  MCRS  is 


{ 


N-  1 
2 


N-  3 
2 


..,  —1,0,1,...., 


N-  3  N  -  1  \ 
2  ’  2  / 


In  what  follows,  we  associate  k  with  the  frequency  slot,  £.  with  the  time 

slot,  and  m  with  code  word  index.  The  following  properties  hold  for  N  prime 

prop.  1)  For  each  code  word  (m  j*  0) ,  each  frequency  slot  and 
each  time  slot  occurs  once  and  only  once.  Thus, 
both  form  CRS's. 

prop.  2)  Each  unequal  pair  of  code  words  (m^n^)  has  one, 
and  only  one,  intersection  Mod  N. 

Property  (2)  is  valid  for  any  horizontal  shift  of  one  or  both  code  words. 

Thus, 


k^(m^,£.)  ■  m^(il-b^)  (Mod  N) 

and 

k2(m2,S,)  **  m2(£,-b2)  (Mod  N)  m^  i  m2 

have  one,  and  only  one,  intersection  for  any  pair  (b^,b2).  This  is  also 
true  for  any  vertical  shifts  since  a  horizontal  shift  of  b  units  is 
equivalent  to  a  vertical  shift  of  mb  units. 
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prop.  3)  For  each  pair  of  unequal  code  words  and  for  any 

horizontal  shift,  the  frequency  and  time  differences 
each  form  a  CRS. 

prop.  4)  For  two  equal  code  words  and  fixed  shift,  the  fre¬ 
quency  and  time  difference  belong  to  the  same  RC. 

The  frequency  and  time  differences  for  all  possible 
shifts  form  a  CRS.  For  a  shift  of  b  time  slots,  all 
frequencies  are  mb(Mod  N)  slot  apart. 

This  completes  the  review  necessary  for  the  sequel. 


IV.  CORRELATION  PROPERTIES 


We  shall  first  energy  normalize  each  code  word.  In  each  segment  of 
length  T/N  seconds,  we  have 

j(V  +  ek) 


fk(t)  =  Ae 


0  <  t  <  T/N 


(3) 


The  energy  of  each  segment  is 

T/N 


fk(t)f*(t)dt  =  a2t/n 


Thus,  since  there  are  N  segments,  we  have 

N 


l 


E,  >AT 
k 


but  E  *  1  thus  A  “  1/ /~T  ".  Thus,  the  rectangular  pulse  has  amplitude  1//  T  , 

Now  we  calculate  the  maximum  crosscorrelation  between  two  segments 

T/N 


mn  T 


j(w  -u  t  j  9  -  6) 

J  m  n  J  m  n  , 
e  e  dt 


-8- 


October  22,  1980 
ELT : cac 


jc% -»„»/» 

1  e  -  1  J  m  n 


U)  -  U) 

m  n 


Taking  absolute  magnitudes  and  observing  that  the  numerator  must  be  no 
larger  than  2,  we  have 


l<*>B-<dnl  *  m  ^  n 


1/N  ,  m  =  n 


Recall  that 


Wm  "  “o  +T  *  ”  =  0,1,. ...N  -  1 


m  n  N 


so  that 


|Cmn|  <_  I  T  |m  -  n|  B 


,  m  +  n 


1/N  ,  m  =  n 


In  order  to  combine  the  results  to  find  the  upper  bound  for  the  cross¬ 
correlation  function,  we  refer  to  Figure  2  in  which  two  code  words  are 
shown  with  an  arbitrary  delay. 


fL  fR 


Figure  2.  Computation  of  Cross- 
Correlation  of  Two 
Code  Words 


•■'y  *  }  '*T  v  -- 
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Suppose  the  lower  pulse  is  a  reference  pulse.  Each  frequency  slot  will, 

in  general,  multiply  two  slots  in  the  upper  pulse.  Thus,  there  will  be 

contributions  from  the  left  portions  (f  's)  and  right  portions  (f's). 

L  K 

We  consider  these  separately.  First  by  property  (2),  if  f,  =  f  ,  we 

^  lj 

have  —  times  the  overlapping  time,  Tq.  Obviously  since  there  will  also 
be  one  other  position  for  which  a  different  slot  will  match  frequencies 
with  right  portions,  and  that  contribution  is  T/N  -  Tq  times  ^  ,  the 
total  contribution  is  For  each  of  the  other  frequencies,  by  property 

(3),  the  difference  will  form  a  CRS.  An  upper  bound  is  clearly  obtained 
if  we  choose  the  MCRS  since  this  is  the  worst  case.  Taking  into 
consideration  left  and  right  segments,  we  have  that 


N  -  1 
2 

,c  I  <  I  +  2  V  -2N 

^cc  -  N  T|kjB 


ki*  0 


N-  1 
2 

|  1  ,  8N  \  1 

lCcJ  -  N  +  BT  k 

k  =  1 


N  -  1 
2 

Thus  '..-M  I  £ 

k  =  1 


(5) 


The  sum  in  Equation  (5)  is  upper  bound  by 


N  -  1 


k  =  1 


r 


In 


N  -  1 
2 


N  >  3 
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so  that  as  a  further  bound  we  have 


c«>  -s  +  f  [1+t*(sr']] 

Figure  3  shows  a  family  of  curves  for  C(N)  parameterized  by  BT.  Although 
these  curves  are  shown  for  all  N,  they  are  only  valid  for  N  prime.  Also, 
observe  that  each  curve  has  a  value  of  N  for  which  it  is  minimum.  In 
order  to  obtain  this  we  differentiate  C(N)  with  respect  to  N,  set  the 
results  equal  to  zero,  and  solve  for  BT  yielding 

<BT>opt  ■ 41i!  [frr + 


N  - 


1 


(7) 


which  is  asymptotically 


(BT)opt~4N2C?  +  £n(N/23 

(N -*■<») 

Substituting  Equation  (7)  into  Equation  (6),  we  obtain 


C(N) 


min 


1 

N 


1  + 


1  +  in 

N-  1 ' 
2 

_1 

2N  -  1  .  „ 
N-l  +X,n 

N- 

2 

V 

1 

rj 

(8) 


(9) 


which  has  an  asymptotic  expression 


C(N) 


min 


1/2-1  < 
8.  n  (N  /  2 )  _  - 


2 

N  ’ 


(N  ■+  °°) 


(10) 


Thus,  we  conclude  that  if  BT  grows  as  Equation  (8),  then  the  uniform  upper 
bound  between  all  elements  of  the  code  set  goes  to  zero;  that  is. 


C(N) 


min 


2 

N  * 
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Figure  3.  Correlation  Bounds  of  the  Code  Words 
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AUTOCORRELATION  PROPERTIES 

Property  four  Indicates  that  after  T/N  seconds  the  frequencies  are 
all  1  slot  (Mod  N)  displaced.  Thus,  we  may  bound  the  autocorrelation  by 
2N  times  the  correlation  for  two  segments,  1  slot  apart.  Thus, 


<  2N  — 
-  _BT_ 


and  if  (BT) 

opt 

(11),  we  have 


is  chosen  from  Equation  (7)  and  substituted  into  Equation 


(12) 


It  is  clear  that  a  trade-off  is  involved.  If  (BT)  „  is  chosen,  the 

opt 

autocorrelation  bound  goes  to  zero  very  slowly,  i.e.,  as  l/£n(N);  whereas, 
if  BT  grows  faster,  say  as  N3,  then  we  obtain  better  sidelobe  properties  of 
the  autocorrelation  but  may  suffer  in  the  crosscorrelation  bound. 


CONCLUSIONS  AND  EXTENSIONS 

We  have  established  that,  for  the  class  of  time  frequency  hop  codes 
defined  here,  we  can  uniformly  upper  bound  the  crosscorrelation  function 
between  any  two  elements  of  the  code  set.  We  have  shown  that  the  bound 
can  be  minimized  by  appropriate  choice  of  BT  product,  which  insures  that 
the  crosscorrelation  function  bound  will  diminish  as  2/N  for  N  large.  It 
should  be  noted  that  the  crosscorrelation  functions  themselves  could  still 
remain  low  since  only  the  bound  could  be  increasing.  However,  several 
computer  simulations  have  convinced  the  author  that  the  bound  is  generally 
tight. 
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In  this  paper,  we  have  results  for  finite  energy  signals  only.  Yet 
the  inclusion  of  finite  average  power  signal,  defined  as  periodic 
extensions  of  the  signals  presented,  is  obvious.  With  the  modifications 


that 


T 

r 


VT)  ‘  T 


f1(t)f2(t  +  x)dt 


o 


and  normalization  constant  in  Equation  (13)  set  to  1,  the  same  bounds  are 
valid  since  the  bounding  procedure  assumed  all  possible  frequency  differences 
were  present  for  all  delays. 
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